function getDenFunValue()
      global dfValue;
      ndist = textread('ndist.txt');
       len = size(ndist,1);
       x = 0:0.01:90;
       border = length(x);
       y = zeros(size(x));
       for i = 1:len
           u = ndist(i,1);
           %sigma = ndist(i,2);
           sigma = 1;
           %sigma = ndist(i,1)/30;
           %sigma = 5;
           for j = 1:border
              y(j) = y(j) +  (1/(sqrt(2*pi)*sigma))*exp(-(x(j)-u)^2/(2*sigma^2)); 
           end
       end
       z = trapz(x,y);
       % display(y);
       dfValue = y/z;
end